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Abstract. Major ingredients of an isospin and momentum dependent transport model for nuclear reactions induced by 
radioactive beams are outlined. Within the IBUU04 version of this model we study several experimental probes of the equation 
of state of neutron-rich matter, especially the density dependence of the nuclear symmetry energy. Comparing with the recent 
experimental data from NSCL/MSU on isospin diffusion, we found a nuclear symmetry energy of E sym (p) ~ 31.6(p/po) 
at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at 
supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made. 



INTRODUCTION 

The Equation of State (EOS) of isospin asymmetric nuclear matter can be written within the well-known parabolic 
approximation, which has been verified by all many-body theories, as 

E(p,8) = E(p,S = 0)+E sym (p)d 2 + 0(S A ), (1) 

where 8 = (p„ — Pp)/{Pp + Pn) is the isospin asymmetry and E sym (p) is the density-dependent nuclear symme- 
try energy. The latter is very important for many interesting astrophysical problems|l], the structure of radioactive 
nucleiJ^EI and heavy-ion reactions0,SS0]- Unfortunately, the density dependence of symmetry energy E sym (p), 
especially at supranormal densities, is still poorly known. Predictions based on various many-body theories diverge 
widely at both low and high densities. In fact, even the sign of the symmetry energy above 3po remains uncertain! 8]. 
Fortunately, heavy-ion reactions, especially those induced by radioactive beams, provide a unique opportunity to pin 
down the density dependence of nuclear symmetry energy in terrestrial laboratories. Significant progress in determin- 
ing the symmetry energy at subnormal densities has been made recently both experimentally and theoretically! 
High energy radioactive beams to be available at GSI and RIA will allow us to determine the symmetry energy at 
supranormal densities. 

To extract information about the EOS of neutron-rich matter from nuclear reactions induced by radioactive beams, 
one needs reliable theoretical tools. Transport models are especially useful for this purpose. Especially for central, en- 
ergetic reactions at RIA and GSI, transport models are the most useful tool for understanding the role of isospin degree 
of freedom in the reaction dynamics and extract information about the EOS of neutron-rich matter. Significant pro- 
gresses have been made recently in improving semi-classical transport models for nuclear reactions. While developing 
practically implementable quantum transport theories is a long term goal, applications of the semi-classical transport 
models have enabled us to learn a great deal of interesting physics from heavy-ion reactions. In the following we out- 
line the major ingredients of an isospin and momentum dependent transport model applicable for heavy-ion reactions 
induced by both stable and radioactive beams 1 11). This model has been found very useful in understanding a number 
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of new phenomena associated with the isospin degree of freedom in heavy-ion reactions. Based on applications of this 
model, we highlight here the most recent progress in determining the symmetry energy at subnormal densities and 
present our predictions on several most sensitive probes of the symmetry energy at supranormal densities. 




IBUU04 VERSION OF THE MODEL 

Crucial to the extraction of critical information about the E sym (p) is to compare experimental data with transport 
model calculations. We outline here the major ingredients of the version IBUU04 of an isospin- and momentum- 
dependent transport model for nuclear reactions induced by radioactive beams 1 11]. The single nucleon potential is one 
of the most important inputs to all transport models. Both the isovector (symmetry potential) and isoscalar parts of 
this potential should be momentum dependent due to the non-locality of strong interactions and the Pauli exchange 
effects in many-fermion systems. In the IBUU04, we use a single nucleon potential derived from the Hartree-Fock 
approximation using a modified Gogny effective interaction (MDI)| 12], i.e., 

U(p,8,p,z,x) = A u {x)^+Ai{x)^ 

Po Po 



po J ' l + (p~p') 2 /A 2 
2C r,t' f ,3 / /V(r,p') 

+ -^rJ d p i+(p- wm- () 

In the above z = 1/2 (—1/2) for neutrons (protons) and t/ t'; (7 = 4/3; fx{f,p) is the phase space distribution 
function at coordinate r and momentum p. The parameters A u (x),Ai(x),B,C> tit ,C T T i and A were obtained by fitting 
the momentum-dependence of the U(p,8,p, T,x) to that predicted by the Gogny Hartree-Fock and/or the Brueckner- 
Hartree-Fock calculations, the saturation properties of symmetric nuclear matter and the symmetry energy of 30 MeV 
at normal nuclear matter density po = 0.16 fm~ 3 | 12]. The incompressibility Kq of symmetric nuclear matter at po is 
set to be 2 1 1 MeV. The parameters A u (x) and A \ (x) depend on the x parameter according to 

2Z? 2Z? 
A u (x) = -95.98 -jc -, A,{x) = -120.57 +x -. (3) 



The parameter x can be adjusted to mimic predictions on the E sym (p ) by microscopic and/or phenomenological many- 
body theories. The last two terms contain the momentum-dependence of the single-particle potential. The momentum 
dependence of the symmetry potential stems from the different interaction strength parameters C TX i and C T T for a 
nucleon of isospin T interacting, respectively, with unlike and like nucleons in the background fields. More specifically, 
we use Cuniike — —103.4 MeV and Q/fe = —11.7 MeV. As an example, shown in Fig. 1 is the density dependence of 
the symmetry energy for x = —2, —1,0 and 1. 




FIGURE 2. Symmetry potential as a function of momentum and density for MDI interaction with x = — 1. 

Systematic analyses of a large number of nucleon-nucleus and (p,n) charge exchange scattering experiments at 
beam energies below about 100 MeV indicate undoubtedly that the symmetry potential at po, i.e., the Lane potential, 
decreases approximately linearly with increasing beam energy according to Uume = a — bE^ where a ~ 22 — 34 
MeV and b ~ 0.1 — 0.2 11131 1 1 411 . This provides a stringent constraint on the symmetry potential. The potential in 
eq|2]meets this requirement very well as seen in Fig. 2 where the symmetry potential (U„ — U p ) /28 as a function of 
momentum and density for the parameter x = — 1 is displayed. 

One characteristic feature of the momentum dependence of the symmetry potential is the different effective masses 
for neutrons and protons in isospin asymmetric nuclear matter, i.e., 



h 2 k dk 



k=k t x 



where k F z is the nucleon Fermi wave number. With the potential in eq. |2j since the momentum-dependent part 
of the nuclear potential is independent of the parameter x, the nucleon effective masses are independent of the x 
parameter too. Shown in Fig. 3 are the nucleon effective masses as a function of density (upper window) and isospin 
asymmetry (lower window). It is seen that the neutron effective mass is higher than the proton effective mass and 
the splitting between them increases with both the density and isospin asymmetry of the medium| 11]. We notice here 
that the momentum dependence of the symmetry potential and the associated splitting of nucleon effective masses in 
isospin asymmetric matter is still highly conto versial ll 1 5l 1 1 fJl . The experimental determination of both the density and 
momentum dependence of the symmetry potential is required. Please see also ref.| 17] on this point. 

Since both the incoming current in the initial state and the level density of the final state in nucleon-nucleon (NN) 
scatterings depend on the effective masses of colliding nucleons in medium, the in-medium nucleon-nucleon cross 
sections are expected to be reduced by a factor 

^medium i a ^ _ [^^j jX NN ) 2 (5) 

where Unn and jX^ N are the reduced mass of the colliding nucleon pairs in free-space and in medium, respecfively fill 
[H ineffective' 

masses is consistent with predictions based on more microscopic many -body theories! 21]. Thus, be- 
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FIGURE 3. Nucleon effective masses in asymmetric matter as a function of density (upper window) and isospin asymmetry 
(lower window). 

cause of the reduced in-medium nucleon effective masses and their dependence on the density and isospin asymmetry 
of the medium, the in-medium NN cross sections are not only reduced compared to their free-space values, the nn and 
pp cross sections are split and the difference between them grows in more asymmetric matter as shown in Fig. [4] The 
in-medium NN cross sections are also independent of the parameter x. The isospin-dependence of the in-medium NN 
cross sections is expected to play an important role in nuclear reactions induced by neutron-rich nuclei|22] 

Other details, such as the initialization of colliding nuclei in phase space, Pauli blocking, etc can be found in our 
earlier publication 1 4ll5l [Till . 



APPLICATION OF THE MODEL 

The model outlined above have many applications in heavy-ion reactions induced by both stable and radioactive 
beams. In this section, we illustrate several examples of studying the role of isospin degree of freedom in the reaction 
dynamics and extracting the density dependence of nuclear symmetry energy. Besides probes of symmetry energy 
at subnormal densities, several probes sensitive to the high density behavior of the symmetry energy have been also 
proposed largely based on transport model simulations. 
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Probing the symmetry energy at subnormal densities with isospin diffusion 

Tsang et al.|9] recently studied the degree of isospin diffusion in the reaction 124 Sn + 112 Sn by measuring|23] 

R _ 2Xi2 *Sn+ il2 Sn ~ Xl24 Sn+ l24 Sn ~ Xn2 Sn+ ll2 Sn ,~ 
X ™Sn+ m S n - Xn2 Sn+ n2 Sn 

where X is the average isospin asymmetry (5) of the 124 Sn-like residue. The data are indicated in Fig. 5 together with 
the IBUU04 predictions about the time evolutions of Rj and the average central densities calculated with x = — 1 using 
both the MDI and the soft Bertsch-Das Gupta-Kruse (SBKD) interactions. It is seen that the isospin diffusion process 
occurs mainly from about 30 fm/c to 80 fm/c corresponding to an average central density from about 1.2po to 0.3po- 
The experimental data from MSU are seen to be reproduced nicely by the MDI interaction with x = — 1, while the 
SBKD interaction with x = — 1 leads to a significantly lower /?,■ value| 10]. 

Effects of the symmetry energy on isospin diffusion were also studied by varying the parameter x\ 10]. Only with 
the parameter x — — 1 the data can be well reproduced. The corresponding symmetry energy can be parameterized as 
E sy m{p) ~ 31.6(p/po) 105 . In the present study on isospin diffusion, only the free-space NN cross sections are used 
and thus effects completely due to the different density dependence of symmetry energy are investigated. As the next 
step we are currently investigating effects of the in-medium NN cross sections on the isospin diffusion. 
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FIGURE 5. The degree of isospin diffusion as a function of time with the MDI and SBKD interactions. The corresponding 
evolutions of central density are also shown. 



Isospin asymmetry of dense matter formed in high energy heavy-ion reactions 

What are the maximum baryon density and isospin asymmetry that can be achieved in central heavy-ion collisions 
at the highest beam energy expected at RIA? This is an interesting question relevant to the study of the EOS of 
asymmetric nuclear matter. To answer this question we show in Fig. 6 the central baryon density (upper window) 
and the average (n/p) p > Po ratio (lower window) of all regions with baryon densities higher than po in the reaction of 
132 5n + 124 Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. It is seen that the maximum 
baryon density is about 2 times normal nuclear matter density. Moreover, the compression is rather insensitive to 
the symmetry energy because the latter is relatively small compared to the EOS of symmetric matter around this 
density. The high density phase lasts for about 15 fm/c from 5 to 20 fm/c for this reaction. It is interesting to see 
that the isospin asymmetry of the high density region is quite sensitive to the symmetry energy. The soft (e.g., x = 1) 
symmetry energy leads to a significantly higher value of (n/p) p > Po than the stiff one (e.g., x = —2). This is consistent 
with the well-known isospin fractionation phenomenon. Because of the E sym (p)8 2 term in the EOS of asymmetric 
nuclear matter, it is energetically more favorable to have a higher isospin asymmetry 8 in the high density region with 
a softer symmetry energy functional E sym (p). In the supranormal density region, as shown in Fig. 1, the symmetry 
energy changes from being soft to stiff when the parameter* varies from 1 to -2. Thus the value of (n/p) p > Po becomes 
lower as the parameter x changes from 1 to -2. It is worth mentioning that the initial value of the quantity (n/p) p > Po 
is about 1.4 which is less than the average n/p ratio of 1.56 of the reaction system. This is because of the neutron- 
skins of the colliding nuclei, especially that of the projectile l32 Sn. In the neutron-rich nuclei, the n/p ratio on the 
low-density surface is much higher than that in their interior. It is clearly seen that the dense region can become either 
neutron-richer or neutron-poorer with respect to the initial state depending on the symmetry energy functional E sym (p ) 
used. 



Pions yields and % /n + ratio as a probe of the symmetry energy at supranormal densities 

At the highest beam energy at RIA, pion production is significant. Pions may thus carry interesting information 
about the EOS of dense neutron-rich matter[24, 25]. Shown in Fig. 7 are the n~ and n + yields as a function of the 
x parameter. It is interesting to see that the n~ multiplicity depends more sensitively on the symmetry energy. The 
n~ multiplicity increases by about 20% while the n + multiplicity remains about the same when the x parameter is 
changed from -2 to 1. The multiplicity of n~ is about 2 to 3 times that of n + . This is because the n~ mesons are 
mostly produced from neutron-neutron collisions. Moreover, with the softer symmetry energy the high density region 
is more neutron-rich due to isospin fractionation! 25]. The K~ mesons are thus more sensitive to the isospin asymmetry 
of the reaction system and the symmetry energy. However, one should notice that it is well known that pion yields are 
also sensitive to the symmetric part of the nuclear EOS. It is thus hard to get reliable information about the symmetry 
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FIGURE 6. Central baryon density (upper window) and isospin asymmetry (lower window) of high density region for the reaction 
of 132 5n + 124 Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 
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FIGURE 7. The it~ and 7t + yields as functions of the x parameter. 




energy from n~ yields alone. Fortunately, the n~ j% + ratio is a better probe since statistically this ratio is only sensitive 
to the difference in the chemical potentials for neutrons and protons Eol . This expectation is well demonstrated in Fig. 
8. It is seen that the pion ratio is quite sensitive to the symmetry energy, especially at low transverse momenta. Thus, 
it is promising that the high density behavior of nuclear symmetry energy E sym (p) can be probed using the n~ /n + 
ratio. 



Isospin fractionation and n-p differential flow at RIA and GSI 

The degree of isospin equilibration or translucency can be measured by the rapidity distribution of nucleon isospin 
asymmetry 8f ree = (N„ —N p )/(N„ +N P ) where N„ (N p ) is the multiplicity of free neutrons (protons)|27]. Although 
it might be difficult to measure directly 8f ree because it requires the detection of neutrons, similar information can be 
extracted from ratios of light clusters, such as, 3 H/ 3 He, as demonstrated recently within a coalescence modell28l l29ll . 
Shown in Fig. 9 are the rapidity distributions of 8f ree with (upper window) and without (lower window) the Coulomb 
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FIGURE 8. The K~ /n + ratio as a function of transverse momentum. 



potential. It is interesting to see that the 8f ree at midrapidity is particularly sensitive to the symmetry energy. As 
the parameter x increases from —2 to 1 the 8f ree at midrapidity decreases by about a factor of 2. Moreover, the 
forward-backward asymmetric rapidity distributions of 8f ree with all four x parameters indicates the apparent nuclear 
translucency during the reaction|30[. 




FIGURE 9. Isospin asymmetry of free nucleons with and without the Coulomb force. 



Another observable that is sensitive to the high density behavior of symmetry energy is the neutron-proton differen- 
tial flowEl 

N(y) 

Fn- P {y)=Win)/N{y), (7) 
1=1 

where Wj = 1(— 1) for neutrons (protons) and N(y) is the total number of free nucleons at rapidity y. The differential 
flow combines constructively effects of the symmetry potential on the isospin fractionation and the collective flow. 
It has the advantage of maximizing the effects of the symmetry potential while minimizing those of the isoscalar 
potential. Shown in Fig. 10 is the n-p differential flow for the reaction of li2 Sn + 124 Sn at a beam energy of 400 
MeV/nucleon and an impact parameter of 5 fm. Effects of the symmetry energy are clearly revealed by changing the 
parameter x. 
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FIGURE 10. Neutron-proton differential flow at RIA and GSI energies 



CONCLUSIONS 

Transport models are powerful tools for investigating especially central heavy-ion reactions induced by neutron-rich 
nuclei. Applications of these models will help us understand the isospin dependence of in-medium nuclear effective 
interactions. Comparing with the experimental data, we can extract the isospin dependence of thermal, mechanical 
and transport properties of asymmetric nuclear matter playing important roles in nuclei, neutron stras and supernovae. 
Currently, the most important issue is the density dependence of the nuclear symmetry energy. The latter is very 
important for both nuclear physics and astrophysics. Significant progress has been made recently by the heavy- 
ion community in determining the density dependence of the nuclear symmetry energy. Based on transport model 
calculations, a number of sensitive probes of the symmetry energy have been found. The momentum dependence 
in both the isoscalar and isovector parts of the nuclear potential was found to play an important role in extracting 
accurately the density dependence of the symmetry energy. Comparing with recent experimental data on isospin 
diffusion from NSCL/MSU, we have extracted a symmetry energy of E sym (p) « 31.6(p /po) 1 ' 05 at subnormal densities. 
It would be interesting to compare this conclusion with those extracted from studying other observables. More 
experimental data including neutrons with neutron-rich beams in a broad energy range are needed. Looking forward 
to experiments at RIA and GSI with high energy radioactive beams, we hope to pin down the symmetry energy 
at supranormal densities in the near future. Theoretically, the development of a practically implementable quantum 
transport theory for nuclear reactions induced by rdioactive beams remains a big challenge. 
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